Electrostatic and electromagnetic instabilities associated with electrostatic shocks: 

two-dimensional particle-in-cell simulation 
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A two-dimensional electromagnetic particle-in-cell simulation with the realistic ion-to-electron 
mass ratio of 1836 is carried out to investigate the electrostatic coUisionless shocks in relatively high- 
speed (~ 3000 km s~^) plasma flows and also the influence of both electrostatic and electromagnetic 
instabilities, which can develop around the shocks, on the shock dynamics. It is shown that the 
electrostatic ion-ion instability can develop in front of the shocks, where the plasma is under counter- 
streaming condition, with highly oblique wave vectors as was shown previously. The electrostatic 
potential generated by the electrostatic ion-ion instability propagating obliquely to the shock surface 
becomes comparable with the shock potential and flnally the shock structure is destroyed. It is also 
shown that in front of the shock the beam-Weibel instability gradually grows as well, consequently 
suggesting that the magnetic field generated by the beam-Weibel instability becomes important in 
long-term evolution of the shock and the Weibel-mediated shock forms long after the electrostatic 
shock vanished. It is also observed that the secondary electrostatic shock forms in the reflected ions 
in front of the primary electrostatic shock. 



I. INTRODUCTION 

CoUisionless shock is one of the most interesting phe- 
nomena in plasma physics; it dissipates the kinetic energy 
of the plasma flow into the thermal energy and the elec- 
tromagnetic energy not by the Coulomb collision but by 
the collective effect associated with the electric and the 
magnetic fields. For example, the universe is filled with 
hot, tenuous coUisionless plasmas and a variety of coUi- 
sionless shocks are produced due to violent phenomena, 
such as supernova explosions. These shocks are believed 
to accelerate charged particles to high energies to gener- 
ate cosmic rays. 

The electrostatic shock [l|, or the ion-acoustic shock 
0, is also one of the coUisionless shocks and it forms 
in unmagnetized coUisionless electron-ion plasmas if the 
Mach number is not so large and the temperature ratio 
of electrons to ions is relatively large Q. The electro- 
static shocks were observed in various experiments with 
double- plasma devices J3-[g| , with Q-machines [3], with 
photo-ionized plasmas and with laser plasmas [3, . 
Recent experiments with intense lasers also showed a pos- 
sible formation of the electrostatic shock at a high shock 
speed of ^ 1000 km ll|. In space, these shocks 



are observed, for example, in the auroral zone of the 
Earth [l3| as well. In the astrophysical context, there 
are no clear observations of the electrostatic shocks so 
far. However, they can also be driven in the universe at 
a wide range of the plasma flow speed. From the the- 
oretical point of view, the electrostatic shock has been 
investigated with hybrid or par ticle-in-cell (PIC) simula- 
tions extensively [e.g.,[3l,[T3l.ll4|. Recent simulations also 
showed a possible formation of very high Mach number 



* E-mail: ka to-t@ile.osaka-u.ac.jp] 



electrostatic shocks [15|. 

The requirement of the relatively high electron-to-ion 
temperature ratio for the electrostatic shocks is also one 
of the conditions for the electrostatic ion-ion instability 
(iGj - fToj . which is also called the ion/ion acoustic insta- 
bility, to develop. Since there exist the reflected ions in 
front of the electrostatic shocks, this instability can grow 
to generate electrostatic waves there. The wave vector 
of the instability is oblique or even almost perpendicular 
to the streaming direction [2l| and therefore multi- 
dimensional simulations are necessary to investigate the 
effect of the instability. Although most of the simula- 
tions were carried out in one dimension and therefore 
could not deal with the instability, Karimabadi, Omidi 
& Quest [2^ carried out a two-dimensional electrostatic 
hybrid simulation and showed that the electrostatic ion- 
ion instability indeed develops in front of the electrostatic 
shock and affects the structure of the shock significantly. 
(It should be noted that recently Ohira & Takahara |23| 
showed that this instability can also develop in the foot 
region of the magnetized coUisionless shocks in supernova 
remnants.) 

In astrophysical plasmas and laboratory plasmas with 
recent laser facilities, flows of coUisionless plasma whose 
velocity is even faster than 1000 km s~^ can be generated. 
In such high-speed flows, the electromagnetic instabili- 
ties, e.g., the Weibel-type instabilities [IjjIlHl) appears 
to be important as well as the electrostatic instabilities 
and can affect the evolution of the electrostatic shocks. 
The simulation by Karimabadi, Omidi & Quest [2^ was 
however an electrostatic one and the effects of such elec- 
tromagnetic instabilities were therefore not included. In 
addition, since the simulation box size of their simulation 
was not large compared with the structures generated by 
the electrostatic instability, larger scale simulations are 
also desirable. 

In this paper, we investigate the electrostatic shocks 
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propagating at a relatively high speed with a two- 
dimensional electromagnetic PIC simulation. In particu- 
lar, we focus on the influence of the electrostatic ion-ion 
instability and the Weibel-type instabilities on the shock 
formation. 



II. LINEAR ANALYSIS 

The electrostatic shocks can form in two counter- 
streaming collisionless plasmas with a relatively high 
electron-to-ion teniperature ratio as shown by the previ- 
ous simulations In such situations, the plasma be- 
comes counter-streaming condition in front of the shock 
and some instabilities (electrostatic and/or electromag- 
netic) can develop there to affect the shock dynam- 
ics. In this section, we present the linear analysis of 
the electrostatic ion-ion instability and the electromag- 
netic beam-Weibel instability in counter-streaming plas- 
mas with wave vectors in arbitrary direction. 

Consider a counter-streaming plasma consisting of 
Maxwellian beams; each beam of the species s has the 
drift velocity of Vg along the x-axis and the thermal veloc- 
ity of Gs- The non-perturbed (zeroth order) distribution 
function for each species is given by 



exp 



(1) 

where Ug is the number density. For this system, the lin- 
ear dispersion relation of the electrostatic mode is given 
by 



2^fcL(i + C^(C)) = o, 



(2) 



where k is the magnitude of the wave vector, 9 is the 
angle between the wave vector and the z-axis, 



1/2 



and 



Cs = Cs('^, k,d) = {uj/k- Vs cos 9) /as 



(3) 



(4) 



The function Z(() called the plasma dispersion function 
[26j is defined by 



— -^-1/2 



-dz 



(5) 



on the upper half of the complex ^-plane (Im(C) > 0) and 
analytically continued into the lower half-plane. Note 
that the solution of the dispersion relation © with 
Im(w) > means an unstable mode. 

Here, we consider a symmetric Maxwellian beam sys- 
tem consisting of two electron beams (s = e+ and e-) 
and two ion beams (s = i-l- and i-) with 71^+ — rie- = 



n,+ = rii- = Uoo, Ve+ = Vi+ = V and K- = Vi- = -V. 
Figure [1] shows the linear growth rate of the electrostatic 
ion- ion instability [l6l - [l9j . F = Im(a;), numerically cal- 
culated for a symmetrical beam system with a relatively 
high electron-to-ion temperature ratio. For later conve- 
nience, we take the following parameters: mi/rric = 1836, 
V = 0.01c ~ 3000 km s-\ = 11.3 keV (oe = 0.21c), 
Ti = 1.25 keV (ai = 1.6 x lO'^c), and thus Te/T; = 9, 
where c is the speed of light, and Tj are the electron 
and the ion temperatures, respectively. The maximum 
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FIG. 1: Linear growth rate of the electrostatic ion-ion insta- 
bility for the symmetrical beam system of electron-ion plasma 
with the drift velocity of 0.01c and the electron-to-ion temper- 
ature ratio of 9. The horizontal and the vertical axes are the 
wave vector in the a;-direction, k^, and that in the {/-direction, 
ky, respectively. The gray-scale shows the linear growth rate 
normalized to the electron plasma frequency. 



growth rate is obtained asF = 2.8xlO 



at k 



where Ac ~ c^pc the electron 



74°, where Wpc = (47rncoe^/mo)-^/^ is 
the electron plasma frequency defined for the mean elec- 
tron number density Uco', therefore the dominant mode 
is directed almost perpendicular to the streaming direc- 
tion. The corresponding typical wavelength is given by 
A = 27r/fc = 6.7xlO-iA. 
skin depth. 

On the other hand, one of the two electromagnetic 
modes whose electric field lies on the plane made by the 
X-axis and the wave vector can be unstable in this system, 
too; this is the beam-Weibel instability (or the filamen- 
tation instability) [25i], which is distinguished from the 
(ordinary) Weibel instability [l^l that grows in a plasma 
with anisotropic temperatures. The dispersion relation 
of the beam-Weibel instability is given by 



-ike) 



a? 



2(K;Sin( 



-il + CsZiCs))-l 



(6) 

Figure [2] shows the linear growth rate of this instability 
for the same system as in Fig. [T] The maximum growth 
rate is obtained as F = 1.5 x lO^'^Wpc a,t k = 0.15 WpeC~^ 
(the typical wavelength of A = 4.2 x lOUo) and 9 = 90°. 
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Although the growth rate of this beam-Weibel instabil- 
ity is smaller by one order of magnitude than that of the 
electrostatic ion-ion instability, the magnetic field gen- 
erated by this instability can be important during the 
long-term evolution. 
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FIG. 2: Linear growth rate of the electromagnetic beam- 
Weibel instability for the same beam system as in Fig. [1] 



III. SIMULATION 



We perform a two-dimensional electromagnetic PIC 
simulation to investigate the electrostatic shocks together 
with the instabilities associated with the shock. The 
simulation code used is a relativistic, electromagnetic, 
particle-in-cell code with two spatial and three velocity 
dimensions, namely 2D3V code, developed based on a 
standard method described in Birdsall & Langdon "27] ; 
this code was also used to demonstrate the formation 
of the Weibel-mediated coUisionless shocks at relativis- 
tic speed in pair plasma f28| and at nonrelativistic speed 
in electron-ion plasma ^29;]. Furthermore, it was used 
to derive the scaling law for the model experiment in 
laboratory to demonstrate the formation of the Weibel- 
mediated coUisionless shocks in laboratory [sO] as well. 
Thus, our code can deal with the electromagnetic modes 
as well as the (oblique) electrostatic modes. The ba- 
sic equations of the simulation are the Maxwell's equa- 
tions and the (relativistic) equation of motion of parti- 
cles. In the following, the simulation plane is regarded 
as the X — y plane and the z-axis is taken perpendicular 
to the plane. We take as the unit of time and the 



electron skin depth Ae 



-1 

CUJ 



pc 



as the unit of length (for 



example, for nco = 1 cm WpJ^ = 1.77 x 10 s and 
Ae = 5.31 X 10^ cr 
fields are E', = B^, 



The units of electric and magnetic 
c(47rnoome)-^/^. 



A. Initial condition 

Here, we take the following parameters, which is sim- 
ilar to those of the one-dimensional simulation shown in 
Fig. 1 (c) of Forslund & Shonk The ratio of the 

ion mass to the electron mass is 1836. The grid size is 
Nx X Ny — 4096 X 1024 and the number of particles is 
^ 54 per cell per species. The physical size of the simu- 
lation box is Lx X Ly = 160Ac x 40Ac and therefore the 
size of a cell is Aa; — Ay ^ 0.039Ac. The bulk velocity of 
the upstream plasma is = 0.01c ^ 3000km s~^. The 
temperatures of the electrons and the ions are To — 11.3 
keV and Ti — 1.25 keV, respectively. Thus, the temper- 
ature ratio is given by Tc/Ti — 9. The corresponding 
thermal velocities are a^. = 0.21c - 6.3 X lO'* km s'^ for 
the electrons and Oi — 1.6 x 10^'^c for the ions. The ion- 
acoustic speed is Cg = (fceTo/wi)^/^ = 3.5 x lO^'^c, where 
/cb is the Boltzmann constant, and the initial bulk Mach 
number is given by Mbuik = V/Cs = 2.9. The Debye 
lengths of the electrons and of the ions are Adc = 0.15Ac 
and Aoi — 0.05Ao, respectively. Initially, the electric and 
magnetic fields are zero over the simulation box. The 
boundary conditions for both the particles and the elec- 
tromagnetic field are periodic in the y-direction. 

In the simulation, a coUisionless shock is driven ac- 
cording to the injection method, or the piston method. 
There are two walls at the left-hand side (smaller x) and 
the right-hand side (larger x) of the simulation box and 
these walls reflect particles specularly. Initially, both the 
electrons and the ions are loaded uniformly in the region 
between the two walls with a bulk velocity of V in the 
-|-a;-direction. The temperatures of the electrons and the 
ions are equal in the upstream. At the early stage of 
the simulation, particles that were located near the right 
wall were reflected by the wall and then interact with 
the incoming particles, i.e., the upstream plasma. This 
interaction causes some instability and eventually leads 
to the formation of a coUisionless shock. Note that the 
frame of the simulation is the downstream rest frame; we 
observe the propagation of the shock from the right to 
the left in the downstream rest frame. 



B. Results 

Figure [3] shows the time evolution of the x-component 
of the electric field, E^, averaged over the y-direction. 
We confirm that the electrostatic shock surely forms at 
around ujpct ^ 1000 and propagates until Wpgi 2000 
at the constant velocity of Vsh.d = 1.2 x 10~'^c measured 
in the downstream frame; the shock velocity measured in 
the upstream frame is given by Vsh = 1-1 x lO^^c and the 
shock Mach number is estimated as Ms — Kh/cs — 3.2. 

After ojpct = 2000, however, the shock slows down and 
finally fades away. Figure |4] shows snapshots of the ion 
number density at ujpct = 1500, 3000, and 5000. It is 
evident that oblique filamentary structures develop in 
front of the shock where the plasma is under the counter- 
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streaming condition due to the existence of the reflected 
ions. We see that with the development of the fila- 
ments, the upstream plasma is strongly fluctuated and 
the shock structure is significantly modified. A sim- 
ilar structure was also observed in the simulation by 
Karimabadi, Omidi & Quest [12]. Figure [5] shows the 
charge density at Wpci = 3000 and its power spectrum, 
which is directly related with the existence of the electro- 
static modes. Comparing the region where the power is 
large in the power spectrum with that in the growth rate 
map obtained by the linear theory shown in Fig. [TJ it is 
clear that these are generated by the electrostatic ion-ion 
instability. Note that, as is shown later, in fact the num- 
ber density of the reflected ions is smaller than that of 
the incoming ions. The ratio is typically 0.6 [see Fig. [8] 
(d)]. However, the growth rate as a function of and 
ky (not shown here) is not significantly different from 
that of the symmetric beam case shown in Fig. [T] the 



maximum growth rate is obtained as F = 2.4 x 

at I A: I — 7.6 WpcC"^ and 9 — 79°, which is only slightly 

smaller than the symmetric case. 

Figure [6] shows the phase-space plots of the ions and 
electrons at u^^t = 1500, 3000, and 5000. The ther- 
mal velocity of the ions in the upstream region increases 
with time in both x and y directions, while the increase 
in the velocity dispersion of the electrons is relatively 
small. Figure [7] shows the time evolution of the velocity 
distributions of the ions within 143 < x/Xc < 145, just in 
front of the shock front. The ions are heated especially in 
the y direction because the electric field generated by the 
electrostatic ion-ion instability is almost directed in the 
y-direction. Thus, the anisotropy in the ion temperature 
increases with time. 

There are two possible causes for the shock decay: one 
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FIG. 3: Time evolution of the x-component of the electric 
field. Ex, averaged over the j/-direction (gray-scale). The hor- 
izontal and vertical axes are the a;-axis and the time, respec- 
tively. The electrostatic shock forms and propagates during 
1000 < ujpct < 2000. However, after that time, the shock 



is the anisotropy in the ion temperature (or the decrease 
in the electron-to-ion temperature ratio), and the other is 
the fluctuation of the upstream plasma; both are caused 
by the electrostatic ion-ion instability. The time evo- 
lution of the temperatures immediately in front of the 
shock are shown in Fig. [5] together with that of the ratio 
of the reflected ion density to the incoming ion density. 
Here, the temperatures are calculated for the particles 
within the region in front of the shock given by 



a;sh(^) - 3Ac < a; < Xsh(i) - lAo 



(7) 



where we take Xsh{t) = -Vsh^it - lOOOWp;,^) + 148. 5Ae 
and Vsh,d = 1.2 X 10~'^c from Fig. [3l For the ions, we dis- 
tinguish the "incoming ions" that have Vx > and the 
"reflected ions" that have Vx < 0, while the electrons are 
regarded as a single population because the mean (bulk) 
velocity of the electrons are negligible compared to their 
thermal velocity. Since the electrostatic shock is formed 




slows down and finally fades away. 



FIG. 4: Number density of the ions at LUpct = (a) 1500, (b) 
3000, and (c) 5000. The horizontal and vertical axes are the 
X-axis and the y-axis in the electron skin depth, respectively. 
The electrostatic ion-ion instability, which generates the fila- 
mentary structure, develops with time. 
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at ojpet ^ 1000, the temperatures before that time are 
meaningless as those in front of the shock. We see that 
the electrons are heated in the x and the y directions in 
front of the shock immediately after the shock formation 
and then each component of the electron temperature 
keeps almost constant. On the other hand, the temper- 
atures in the y-direction for both the incoming ions and 




the reflected ions rapidly increase with time within the 
period 1800 < LUp^t < 2300, which coincides with the 
period of the shock decay, while the other components 
do not change significantly. The ion heating in the y- 
direction would be caused by the electrostatic ion-ion in- 
stability and results in the anisotropic temperatures in 
the ions as well as the reduction of the electron-to-ion 
temperature ratio. 

To clarify whether the anisotropic temperature in 
the ions is the cause of the shock decay, we carried 
out another (quasi) one-dimensional simulation with an 
anisotropic ion temperature. The simulation code used is 
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FIG. 7: Velocity distributions of the ions in the v^-Vy space 
within 143 < x/Xc < 145 at Wpet = (a) 1500, (b) 3000, and 
(c) 5000. The heating in the j/-direction is remarkable. 



FIG. 5: (a) The charge density at Wpci = 3000 and (b) its 
power spectrum. The horizontal and vertical axes in (b) are 
kx and ky in cjpcC"'^, respectively. The portion where the 
electrostatic mode develops shown in the panel (b) is well 
agree with those obtained from the linear analysis shown in 

Fig.m 
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FIG. 6: Phase-space plots of the ions and electrons at uupat = 
(a) 1500, (b) 3000, and (c) 5000: from top to bottom, the x-v^ 
and x-Vy distributions of the ions and those of the electrons 
are shown. In the upstream of the shock, the ion heating is 
evident while the electron heating is small. 




FIG. 8: Time development of the temperature in each direc- 
tion immediately in front of the electrostatic shock normalized 
to the initial ion temperature Tio: (a) the electrons, (b) the 
incoming ions, and (c) the reflected ions. The solid curves, 
the dashed curves, and the dotted curves show the x, y, and 
z components, respectively. It is remarkable that the y com- 
ponents of the temperatures of both ion populations rapidly 
increase within the period 1800 < ujpct < 2300. The panel 
(d) shows the time development of the number density ratio 
of the reflected ions to the incoming ions, ni-/ni+. 
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the same two-dimensional code but we take only 8 grids 
in the y-direction so that the simulation is essentially one- 
dimensional; in such a simulation, since the electrostatic 
ion-ion instability does not develop, we can investigate 
the pure effect of the anisotropy in the ion temperature 
on the shock formation. For the shock formation, the 
condition of the incoming ions is more important than 
that of the reflected ions and Fig.|5](b) shows the typical 
anisotropic temperature of the incoming ions is given by 
Ti^y ~ 2Ti_2. around the time of the shock decay. Here, 
as a more severe case, we take T^ y — 4Ti x keeping the 
other parameters unchanged as in the two-dimensional 
simulation, that is, Tc/Ti^x — 9. The time evolution of 
the electric field Ex is shown in Fig. [9] (a) . We found that 
the electrostatic shock forms under even this condition, 
where the shock speed measured in the upstream frame 
is Vsh ~ 1.2 X lO^^c and the Mach number is given by 
Ms - 3.5. 

To confirm that the decrease in the electron-to-ion tem- 
perature ratio with isotropic ion temperature is also not 
the cause of the shock decay, we performed another quasi 
one-dimensional simulation with an isotropic ion temper- 
ature but the lower electron-to-ion temperature ratio of 
Te/Ti = 4. The result is shown in Fig. [S] (b) showing 
that the shock is also formed with even this temperature 
ratio at almost the same shock velocity and the Mach 
number as the anisotropic case. Thus, the anisotropic 
temperature in the ions nor the decrease in the tempera- 
ture ratio are not the cause of the shock decay observed 
in the two-dimensional simulation. 





FIG. 9: Time evolutions of the a;-component of the electric 
field. Ex, as in Fig. [3] obtained from quasi one-dimensional 
simulations (see text): (a) for the anisotropic ion temperature 
Ti^y/Ti^x = 4 and TejTi^^ — 9 and (b) for the isotropic ion 
temperature with Te/T\ = 4. The electrostatic shocks form 
and propagate without decaying in both cases. 

For completeness, we also performed another quasi 
one-dimensional simulation with the same parameters as 
the two-dimensional one, that is, T^/Ti = 9, except the 
grid number in the y-direction. We found that the elec- 
trostatic shock forms and propagates without decaying in 
this case as well (not shown) confirming that the shock 
decay is a multi-dimensional effect. 

Figure [TU] shows the electrostatic potentials, the num- 



ber density of the incoming ions, and that of the reflected 
ions at ujpct = 2000 for the two-dimensional simulation. 
The potentials are normalized by the upstream ion ki- 
netic energy {miV'^/2). It is evident that the large po- 
tential fluctuation, which reaches even a half of the shock 
potential, exists in front of the shock due to the electro- 
static ion-ion instability. Because of this electrostatic po- 
tential, the densities of both incoming and reflected ions 
are significantly disturbed where the fluctuation pattern 
almost traces the waves of the electrostatic ion-ion in- 
stability. In particular, it can be seen that the portions 
where the incoming ions enter the "shock front" are dif- 
ferent from those where the ions are reflected toward up- 
stream from the "front". Figure [TT] shows the y-Vy phase 




FIG. 10: Close-ups of (a) The electrostatic potentials, (b) the 
number density of the incoming ions (ni+), and (c) that of the 
reflected ions (ni_) at ujpct = 2000 for the two-dimensional 
simulation. In (a), the electrostatic potential in front of the 
shock, which is generated by the electrostatic ion-ion insta- 
bility, reaches even about a half of the shock potential. 

space distribution of the ions around x 145. 5Ac. The 
feature of the ion trapping [l^, [lH is clear and there- 
fore the fluctuation of the ion densities is a result of the 
nonlinear evolution of the electrostatic ion-ion instability. 
Thus, it can be concluded that the cause of the decay of 
the electrostatic shock is the large fluctuation in the ion 
density due to the electrostatic ion-ion instability. 

In Fig. [To] (c) , we can also find another interesting fea- 
ture in the refiected ion density that is seen as a dis- 
continuity around x 140Ao, which could also be found 
in Fig. m (a). The x-Vx phase space distribution of the 
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ions is plotted in Fig. [121 As is shown in the panel (a), 
there is a structure in the phase space of the reflected 
ions around x ^ 140 Ac. The panel (b) shows the close- 
up of this structure. We see that some of the reflected 
ions are reflected again there in the opposite direction. 
This structure would also be another electrostatic shock 
in which the upstream region is the larger x side. This 
'secondary' electrostatic shock would be formed as a re- 
sult of that the faster reflected ions run into the slower 
ones from behind. 

As was shown in Sec. HIl the beam-Weibel instabil- 
ity can develop in counter-streaming plasmas with the 
wave vector in the y-direction as well as the electro- 
static ion-ion instability. However, the large anisotropy 
in the ion temperature {Ti^y/Ti^x 2 for the incoming 
ions and Ti^y/T{^x ~ 6 for the reflected ions) caused by 
the electrostatic instability can affect the growth of the 
beam-Weibel instability. In addition, this temperature 
anisotropy itself may also cause the (ordinary) Weibel 
instability that has the wave vector in the x-direction. 
Here, let us consider the linear dispersion relation of 
these two Weibel-type instabilities in the anisotropic 
Maxwellian beam system given by the following zeroth- 
order distribution functions: 



fQ\vX,Vy,V,) 



T-3/2, 



■ exp 



'sx^sy^sz 



{Vx - Vs)' 



y 

sy 



,(8) 




FIG. 11: The y-Vy phase space distribution of the ions within 
145.5 < x/X^ < 145.7 at Up^t = 2000 (c.f. Fig. [TO)). The 
feature of the ion trapping due to the nonlinear evolution of 
the electrostatic ion-ion instability is clear. This trapping 
results in the strong fluctuation in the incoming ions as is 
shown in Fig. llOf bl to destroy the shock structure. 
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FIG. 12: (a) The s-ii^ phase space distribution of the ions 
within < y/Ae < 0.4 at Wpet = 2000. (b) The close-up 
of the distribution of the reflected ions around x ^ 140 Ac. 
The "secondary" electrostatic shock exists in the reflected ion 
beam around x = 140Ae • 



where Ugx, OLsy and Cgz are the anisotropic thermal veloc- 
ities for the respective directions of the species s. In 
the following, we consider the beams are symmetrical 
for simplicity, that is, ne+ = n^- = rii^ — rti_ — rico, 
K+ ^Vi+^V and Vc- ^Vi- = -V, as in Sec. M Un- 
der this condition, the linear dispersion relation of the 
beam-Weibel instability with the wave vector in the y- 
direction is given by 



c.2-(fcc)^+5:< 



2V' 



(i + C^(Cs))-i 



= 0, 



(9) 



where C.s = uj/kagy. On the other hand, the linear dis- 
persion relation of the (ordinary) Weibel instability with 
the wave vector in the a;-direction is given by 



aix 



0, (10) 



where rj^ = (w/fc — Vs)lasx- (Note that when ^ = 0, this 
relation is reduced to that derived by Weibel [IJ] with- 
out the background magnetic field.) Figure [13] shows 
the growth rates of the two modes. We take here 
T\,ylT\,x = 4, T^,z = T\,x and T^Ti^x = 9 as typical val- 
ues and the other parameters are the same as those of the 
initial condition of the two-dimensional simulation. As 
a reference, the growth rate for the isotropic beam case 
with Tc/Ti = 9 discussed in Sec. [iTj is also shown in the 
figure. We see that the growth rate of the beam-Weibel 
instability is reduced compared with that of the isotropic 
beam case due to the temperature anisotropy in the ions. 
Under this condition, the Weibel mode with the wave vec- 
tor in the cc-direction is stable and does not grow. (Of 
course, in the limit — > 0, this mode should be unstable. 
We confirmed that it is unstable for V < 1.15 x 10~^c for 
the same parameters.) 

Figure [Til shows the z-component of the magnetic field, 
Bz, obtained from the two-dimensional simulation at 
ujpct — 1500, 3000, and 5000. We see that the magnetic 
field grows with time with a wavelength comparable with 
that derived by the linear analysis 40Ao), although 
it may be affected by the periodic boundary condition 
because the wavelength is comparable with the system 
length in the y-direction. At this time (wpoi ~ 5000), the 
typical magnetic field strength is ~ 2 x 10~^ and 

the corresponding ion gyro-radius is given by ~ 9 x 10^ Ac; 
therefore, the ions are hardly deflected by this magnetic 
field. However, in the Weibel-type instability, the mag- 
netic field is generated by the current filaments and the 
current filaments coalesce each other to grow larger fila- 
ments with stronger magnetic fields until saturation [see 
e.g.,'32'1. Although, because of the limitation of the sim- 
ulation box, we cannot study the further evolution of the 
filaments in this simulation, in the real world or in much 
larger simulations, the magnetic field would grow much 
stronger with the coalescence of the filaments and, long 
after the electrostatic shock vanished, it would become 
strong enough to form the "Weibel-mediated" shocks, a 
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kind of collisionless shocks that dissipates the upstream 
particle kinetic energy via the magnetic field generated 
by the instability in the shock transition region. 

As is reported in our previous paper of Ref. , we 
have performed a PIC simulation of the Weibel mediated 
shock with the larger flow velocity {V — 0.45c) and the 
smaller mass ratio (mi/mo — 20) with the equal tempera- 
tures for the electrons and the ions. (In that simulation, 
the electrostatic shock and the electrostatic ion-ion in- 
stability do not develop.) To illustrate the formation of 
the Weibel-mediated shock qualitatively, the time evo- 
lutions of the ion number density and the magnetic field 
(z-component) are shown in Fig. [151 The figures show the 
ion current filaments develop with time due to the beam- 
Weibel instability and finally the Weibel-mediated shock 
forms at ujpct = 2000. Note that due to the mass ratio 
and the fiow velocity, the temporal and spatial scales are 
different from those in the system with the real mass ra- 
tio of rrii/mc — 1836 and the slower fiow speed ~ 3000 
km s^^. Returning to the situation of the present study, 
the ion current filaments developing around the electro- 
static shock (Fig. [Ti| would also grow further and finally 
will evolve into the Weibel-mediated shock in the real 
world or in the simulations with much larger spatial and 
temporal scales. 



IV. CONCLUSION 

We have carried out a two-dimensional electromagnetic 
PIC simulation for the case of counter-streaming plasmas 
at a relatively high fiow velocity with a large electron-to- 
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FIG. 13: The linear growth rates of the two Weibel- 
type instabilities in the anisotropic Maxwellian beam system 
(Ti,y/Ti,^ = 4 and Te/Ti,^ = 9): the beam-Weibel instabil- 
ity with the wave vector in the t/-direction (solid curve) and 
the (ordinary) Weibel instability with the wave vector in the 
a;-direction (dashed curve). The dot-dashed curve shows the 
growth rate of the beam-Weibel instability in the isotropic 
Maxwellian beam system as in Sec. |ll] (i.e., Ti^y/Ti^x = 1 
and Tc/Ti = 9) as a reference. The growth rate of the 
beam-Weibel instability is reduced due to the temperature 
anisotropy in the ions. The Weibel instability in the x- 
direction is stable in this case. 



ion temperature ratio of 9. At first we have confirmed 
that the electrostatic shock forms in the early time evo- 
lution as was shown in the previous works [3, [l^ . We 
also confirmed that the electrostatic ion-ion instability 
develops in front of the shock due to the counter-streams 
of the ions 22]. Then, we found that the electric field 
generated by this instability results in the strong fiuc- 
tuation in the ion density and finally leads to destroy 
the shock itself. It was also found that the electromag- 
netic beam-Weibel instability develops much slower than 
the electrostatic instability but it becomes predominant 
in the later time. This suggests the possibility that the 
Weibel-mediated collisionless shock is formed according 
to the scenario shown in Kato & Takabe [l^ long after 
the electrostatic shock disappears. It was also observed 
that the secondary electrostatic shock forms in the re- 
fiected ions in front of the primary electrostatic shock. 
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FIG. 14: The z-component of the normalized magnetic field, 
B^/B,, at cjpct = (a) 1500, (b) 3000, and (c) 5000 for the 
two-dimensional simulation, where -B* = c(47moo'Tio)^''^. The 
magnetic field is generated by the beam-Weibel instability in 
front of the shock and it grows with time during the simula- 
tion. 
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FIG. 15: PIC simulation of the Weibcl-mcdiatcd shock for 
showing the formation of the shock qualitatively with param- 
eters of rrii/me = 20 and V = 0.45c. The left panel shows the 
ion number density normalized by the upstream value, rii/riao, 
and the right panel shows the ^-component of the normalized 
magnetic field, Bz/B,, from top to bottom at iOpct — 500, 
1000, and 2000, respectively. Due to the magnetic field gen- 
erated by the beam-Weibel instability, a kind of coUisionless 
shock, which is called the Weibel-mediated shock, is formed 
at Wpet ~ 2000. 
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